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1 Introduction 

In applied sciences, many model problems take the formulation of fourth order elliptic perturbation 
problems, such as, e.g., the linearized Cahn-Hilliard equation [9l l3ni39ll43ll5TW53] . and the strain gradient 
problem [ininimEraH]. In order for the robust discretisation of such problems, numerical schemes 
that work for both fourth and second order problems are needed. The rectangular Morley element (RM 
element for short in the sequel) scheme is one that falls into this category. The RM element is introduced 
by Wang, Shi and Xu ( [40], 2007) originally for fourth order problem in arbitrary dimension; its a priori 
( [40]). a posteriori ( i) and superconvergence ( |T8]) analysis have been established already. Utilised 
for second order problems, the RM element scheme is of nonconforming type, and its validity for second 
order problems has been pointed out in Shi-Wang m, however, without a formal statement and technical 
proof. In this paper, we will present a technical and complete analysis of the RM element scheme utilised 
on the second order problem in arbitrary dimension. Specifically, for the energy norm of the error, beside 
the standard analysis which lead to a convergence of the scheme which is of 0(h) order on general shape 
regular triangulations, a more careful analysis is given; namely, when the triangulation is divisionally 
uniform (see its precise description in Section[3l), the convergence rate can be of 0(h 15 ) order, and when 
the triangulation is uniform, the convergence rate is 0(h 2 ). For the L 2 norm, an 0(h 2 ) convergence rate 
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can be obtained for convex domain; and in general, this estimate can not be improved. The discussion 
on these uniform triangulations illustrates the convergence analysis being a sharp one. 

As the RM element is of nonconforming type, revealed by the Strang lemma, work has to be spent on 
the analysis of the consistency error. It is well known that once the zero-th order or first order moment 
of the finite element function is continuous across the internal interfaces of the grids, the consistency 
error of first order (see, e.g., [33tIll[2Dl[MlESlSl]) or of second order (see, e.g., [221111]) can be proved. 
Besides, based on the symmetry of the rectangular cells, another standard way can be to constructing 
an internal orthogonal space on every cell, which can lead to a consistency error of first order (see, e.g., 
the Wilson element) or second order (see, e.g., [2HE2] ) by the internal-eliminating technique. However, 
neither of the techniques above works for the RM element directly. Actually, the average of the RM 
element function is not continuous across the internal interfaces, and a direct utilization of the internal- 
eliminating technique for the RM element can lead to a consistency error of first order only, which can 
not explain the high accuracy of the scheme on uniform grids. By the aid of stable decomposition, we 
can indeed divide the consistency analysis on the whole finite element space to an equivalent system of 
subproblems on a big subspace associated with vertex degree of freedoms (DOFs) and a series of smaller 
subspaces each associated with a face DOF. We can thus implement the internal-eliminating technique 
with respect to cells for the big and small subspaces, and generalize the internal-eliminating technique 
onto patches for every small subspace, and arrive at the convergence rate in energy norm with respect 
to various triangulations, namely 0(h ) on general grids and 0(h 2 ) for uniform grids. The application of 
duality argument seems standard, and an error estimate of 0(h 2 ) is achieved in L 2 norm. However, we 
note that people can not find a nontrivial conforming subspace in the RM element space; the bilinear 
clement space, e.g., is not contained in. This implies that we can not expect the L 2 norm of the error 
being one order higher than the energy norm of the error. Indeed, we prove rigorously for uniform 
triangulations that the L 2 norm of error can not be nontrivially higher than 0(h 2 ). All the analysis 
above is carried out in a unified form for arbitrary dimensions, and numerical experiments confirm the 
theoretical results. 

We remark that, a key fact of our analysis is higher consistency accuracy can be expected on uniform 
triangulations. This fact is studied in the context of superconvergence. Also, this higher accuracy analysis 
has been a basis of further study of the patch recovery technique and a posteriori error analysis. We 
refer to [7ll8l ll2lil41l231l25ll36l[55H59l for related discussions. By the general theory of fast auxiliary space 
preconditioning (FASP) ( c.f. [13114611471154] 1. the stable decomposition will also play a fundamental role 
in designing optimal preconditions in future. 

In the sequel, we will use the following standard notation. We use H for a general bounded polyhedral 
domain in K. d (d ^ 2), dtt the boundary of Q, and n = ^ni,n 2 , • • -,nfj the unit outer normal to dCt. 
For a nonnegative integer s, we shall use the usual Sobolev spaces such as H S (K) with the corresponding 
seminorm and norm denoted by | •| S) if and 11 ■ | \ s ,k, respectively. (•, -)k denotes the inner product of L 2 (K). 
When K = f2, we just write |-| s , ||-|| s and (•, •). Given a multi-index a = (ctq, •••, a^), set |a| = X)f=i a i and 
x a = • • -Xj d , Mx £ R d . For a subset B C and a nonnegative integer r, Let P r (B) and Q r (B) be the 

spaces of polynomials on B defined by P r {B) = span{x"||a| < r}, Q r (B) = span{ 2 Q |aq < r, 1 < * ^ d}. 
In this paper, we use C to denote a generic positive constant which may be different at different places. 
Also, following [45], <, > and = denote ^ and = up to a constant, respectively. The hidden 
constants depend on the domain, and, they also depend on the shape-regularity of the triangulation 
when it is involved, but they do not depend on h or any other mesh parameter. 

The remaining of the paper is organized as follows. In Section [2] we present some preliminaries of the 
RM element in any dimension. In Section [3] we study the discretisation scheme of the model problem, 
and we construct the error estimates in energy norm and L 2 norm on general d-rectangular triangulations. 
In Section [4] some numerical examples are presented to demonstrate our theoretical results. Finally, in 
Section [5] some conclusions are given. 
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2 Preliminaries 


2.1 The d-rectangular Morley element 

Let AT C be a d-rectangle, x c = (xi )C , • • • , Xi lC , • • • , Xd, c ) T 6 be the barycenter of fb, and hi the 
half length of K in x, direction, i = 1,2 ,d. Then the d-rectangle can be denoted by 

K — — (xi, * • • , Xi , • * • , Xd ) | Xi — Xi^ c p 1 b SZ d}, (2.1) 

and the vertices a*, 1 ^ i ^ 2 d , of AT are denoted by 

— (* 1,0 ~b ^ 1 ^ 1 ) ’ ' ' ; *^i,c “t“ ^jdhji * * * )Xd,c “t“ ( \idhd ) ; I — 1 ) 1 ^ J ^ d, 1 $C ?^2 . 
Moreover, let and F 2 j (1 < j ^ d) denote the two (d — l)-faces of A' perpendicular to a^- axe as 

F 2 j — 1 — — (^ 1 , ) ) * 5 *^d) | — *^i,C + l^Z^d, f 7 ^ J, — 1 }, 

and F 2 j — {ir (aii, * * * * * * i^Cd) \ Xi — Xi^ c P £,ihi, 1 ^ ^ 1 , 1 ^ i ^ d, f y- j • — 1 }. 




Figure 1 Degrees of freedom of the rectangular Morley element in two and three dimensions. 


The d-rectangular Morley element ( [2333) is defined by the triple ( K,Pm(K),D ), where 

• the geometric shape AT is a d-rectangle; 

• the shape function space is 

Pm(K) := Qi(AT) +span{xf,xf 11 < i ^ d}; (2.2) 


the vector D{v) of degrees of freedom is, for any v £ C 1 (AT) (see Figure [l]), 

£> ( V ) ; = ( V ( aO, ^/ F ^T rfs )- * = 1, • • •, 2 d , j = 1,..., 2d, 


(2.3) 


where n p j is the unit normal vector of (d— l)-face Fj , and \Fj\ denotes measure of (d— l)-face Fj. 


The triple is Pif-unisolvent. Indeed, define 

d d 

t , 

Pi 




i=i 


13 h 
r h 


j x 1 -x 2 , {{ x_ l -x ljl)2 _ m 


j =1 


h'k/'Ek %k,c . Z'k.c 

Q2k—1 = -r( - 7 - + 1 ) (- 7 -- - 1 ), 

4 hk fik 


hk / %k *£fc,< 

“ = “T ( —ST” 


l)( Xk , gfc - c -l) 2 , 

hk 


1 ^ fc < d, 
1 ^ k < d. 


(2.4) 
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Then it can be verified that, with 8ij the Kronecker symbol, 


Pi (pj ) — A 


l<bj<2 a 


Tpn \ T^ ds = °- 1 < i < 2 d , 1 < j < 2d; 


j\ j Fj dn Fj 

Qi( a j) = 0 , 

1 f dqj 

\ F j\ Jf , dn Fj 


1 sj * s; 2d, 1 < j ^ 2 d , 
ds = Sij , 1 ^ i, j ^ 2d. 


(2.5) 


The corresponding interpolation operator fix is then given by 


2 “ 2d x 

n KV = J2Pi v ( a i) + J2 q jTjT\ / 

2=1 j'=l 1"^' •'■Fj 


<9u 

<9n Fi 


ds,Wv £ C^if). 


Denote P$(if) = span{pi, 1 < i ^ 2 d }, P^(if) = span{g i; 1 ^ < 2d}, and Q k M {K) = span{g 2 fc-i, q 2 k}, 1 < 
k < d. Besides, for / a (d — l)-face of if, denote 

pIj(K) := {u G Pm (if) : u vanishes on vertices of if, 


f 

If dn f' 


ds vanishes on any face f C <9if other than /}. 


Further, given / a (d— l)-face, we can construct two d -rectangles K F and K F that share / as a common 
face. Denote w/ := if/, U if_R as the patch associated with /. |df/,| and \Kr\ denote the measure of 
elements if/, and Kr , respectively. If |if/,| = |if//|, then the patch is called uniform. Define 

Pm(^/) := (vfc € L 2 (D) : v/J/r G Pj(j(if), if £ (if/., if/?}, [ ds is continuous on /}. 

Jf dn f 

Evidently, dim(P/ f (A")) = 1. 

2.2 Structural properties of the shape function space 


2.2.1 Local orthogonality of the finite element function 

As the foundation of the theoretical analysis, some facts on local orthogonality or near orthogonality 
have to be revealed. By direct calculation, we obtain Lemmas 12.11 and 12.21 below. On every element if, 
let be the nodal interpolation operator associated with Qi(if) element. 

Lemma 2.1. Let if be a d-rectangle, n A be the outside normal direction of dK and nf be the i-th 
component of n A , i = 1 ,,d. Then 

1. it holds for <j> £ P^(K) and pf £ Pi (if) that 

f Pi{<t> — Hxfynfds = 0, 1 ^ i ^ d; (2.6) 

J dK 


2. it holds for ifj £ Pfj(K) that 


ipnf ds = 0, 1 ^ ^ d. 


I OK 


(2.7) 

Zi-Zi.c < l < d _ 


Proof. Let if be denoted as (12.11) . In order to simplify the presentation, denote = — h 

1. For any <j> £ Pj^(K), we have (f — lP K (f £ span{£j(£ 2 — 1), 1 < j < d}. It can be verified that 

[ &•£, (Of - l)ds = f fkfjif] - l)ds = 0, 1 < i ^ d, 0 < k ± j < d, 

JF2i — \ J i*2i 
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where £ 0 := 1. Now given pf = c 0 + ci£i + • • • + Cd£,d , thus Pi\F 2i _ 1 = <H — c,£j + Efc=o c fc6= an d 
pf |F at = -Q - + Eto c k£ki we have, for i = 1,..., d, 


[ Pitjitj- l)n?ds= [ Pi — l)ds — f Pi^(^-l)ds 
JdK J F 2i - 1 J F 2i 

= [ ( c j€j - f -'6 - 1 )ds - / (Cjtj ~ - 1 )ds 

J Fii —i J F 2i 

[ ( c jtj - c j£i)Zj(Zj -l)ds- f (c^ -Cj&tjit? - l)ds = 0, if * = J, 

J F 2 i —i J F 2 i 


/ c J^(^- 1 ) ds -/ c i£j(£i “ !) ds = °> 

' F 2 i—\ J F 2 i 


if * ^ j. 


( 2 . 8 ) 


2. Let g € Pf^(K), then gl^-i = <?|F 2 i as two functions of (xi, ..., Xj_i, Xj+i, ..., x^), * = 1,..., d. 
Therefore, for i = 1,... ,d, 


/ qnfds = / qds — / gds = 0. 

IdK ■/ F 2 i-\ JF 2i 


(2.9) 


The proof is finished. □ 

Lemma 2.2. Let / be a (d — l)-face, w/ be the patch associated with /, n^f be the outside normal 
direction of dcjf, and be the i-th component of n“t. If is uniform, then it holds for ip £ Pf (ujf), 
Pi f £ Pi(ujf) that 

f p“ f ip'nf f ds = 0, i = l,...,d. (2-10) 

J dco f 

Proof. Let the barycenter of {d — l)-face / be (/i, c; ■ • • , fk.c, ■ ■ • , fd,c) and the half length of Kl and Kr 

in Xi direction be hi , i = 1, 2,..., d. Then, with Xj being orthogonal to /, we can denote 

f — {x — (xi, ‘ " * , Xfc, • • ■ , Xd) | Xfc — f k.c P ^khk? 1 ^ £k ^ 1? 1 P k P d : h f 1 j. (y — 0}, 

P L — {x — (Xl i 5 Xfc 7 i Xd ) | Xk — fk.c P f,k hki 1 ^ 1 5 k P j. 2 ^ ^ 0}, 

and 


P R — {x — (Xl, ■ ■ ■ , Xfz , , Xd) | Xk — fk.cP dpkhk ■ 1 i 1 7 ^ J, 9 ^ f(y -^7 2 } . 

Now, without loss of generality, let = cq P C 1 X 1 + • • • + c^x^ on ujf and ip £ Ph(u)f), namely, 


ip = 


(fj,c hj) + 2 Xj {fj.c hj) _ 

4 hj hj 

~ Vi* + h i) + 1)( ^ ~ + _ 1)2 quKr: 

4 hi hi 


with some atl. Elementary calculus leads to that 


>dK L 


Pi f ipn^ f ds = 


0, 


if j = i, 


— ^Cihja\KR\. otherwise; 


and, 


ldK R 


p 1 f ipn i ’ ds 


0, if j = i, 

^Cihja\Kji\, otherwise. 


Since p^ f and ip are both continuous across /, w/ is uniform and thus \Kl\ = \Kr\, we have 

[ Pi f ipn“ f ds= f Pi f ipn^ f dsp f p^ f ipn^ f ds = 0, i = l,...,d. 
Jdojf J8 Kl JdKR 


□ 


This finishes the proof. 

Lemma 2.3. There exists a constant 6 £ (0, ^) depending on d only, such that, for any K a d-rectangle, 
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1. it holds for nodal basis function q z . q 3 (1 ^ i ^ j ^ 2 d), that 

I(V 9i , Vfc-)*I < e(\\v qi \\l K + HV^-IIiU); ( 2 . 11 ) 

2. it holds for any </> G Pfj(K) and if G Pfj(K) that 

l(V0, W>kl < 0(l|V^||^ + IIV^II^). (2.12) 

Proof. Firstly, direct calculation leads to that, with i ^ j. 

|(Vg», Vqj)ic\ j gi if {i,j} = {2/c — 1, 2k} for some k, 1 < k < d; 

HVftll^ + IIV^II^ [ 0, otherwise. 

Thus (12.111) is proved. 

Secondly, according to the definition, Pf r (K) can be decomposed as 

Pm(K) = Qm(K) 0 • • • © Qm{K). (2.13) 

Moreover, the decomposition is orthogonal with respect to the inner product (V-, V-)k- Actually, dtq = 0 
for q G Q j m (K) with i ^ j. Further, for every k, 1 < k ^ d, Q k M {K) can be decomposed as 

QUK) = span{(j 2 fc-i,i? 2 fe} = Span{q 2 fe -1 + q 2 k} © span{g 2 fe-i - q 2 k}- (2.14) 

Where (g 2 fc -1 +<Z 2 fc) is perpendicular to P^(K), and this decomposition (12.141) is orthogonal with respect 
to both the inner products (d k -,d k ■)k and (V-, V-)k- 

Meanwhile, for every k, 1 ^ k ^ d, we make an orthogonal decomposition of P^(K) as 

Pm(K) = Wk(K) 0 Y k (K), (2.15) 

where W k (K) is perpendicular to (<Z 2 fc-i ~ Q 2 k), and Y k (K) is the orthogonal complementary of W k {K) 
with respect to the inner product {d k qd k -)x- Note that dim(span{(g 2 fe-i — Q'ik)}) = 1 and ((? 2 fc-i — <? 2 fe) 
is not orthogonal to Pfj(K), we have obviously dim(Yfc(A')) = 1. Denote by y k the unique (up to a 
constant) basis function of Y k (K). 

Now, given (j> G P^(K) and if G Pfj(I\), they can be decomposed as 

<t> = <t>'k + <t>'L $k G span{y/j} and <f" K G W k (K), (2.16) 

and 

d d 

= +^k): (2-!7) 

k =1 k=1 

where if k G Q k M {K), if k G span{g 2 fc-i - <? 2 fc} and if k G span{<? 2 fc-i + <? 2 fe}- Then we have 

d d d d 

(v& vvo* = 53(v^,v^ = J2( dk ^ dk ^« = • 

k= 1 k= 1 k= 1 k= 1 

Note that d k <f' k ^ d k if' k unless both of them are zero; actually, the face average of d k y k vanishes for every 
face, and the face average of d k (q 2 k -1 — <? 2 fc) does not vanish for F 2 k -i or We have 

(dk<t>k,d k il>k) € 0k(\\dk<l>k\\l tK + 

with 9 k < 1/2 uniform on span{j/fc} x span{q 2 fc-i — q 2 k}- Therefore, 
d d 

I(V</,WOid = \Y,WkMk)K\ < E>(ll^*llo,* + ll^fcllo,if) 

/e=l fc=l 

d d 

(E Iloilo,A- + E PfcV’fcllo.ic) = 0(||V^||g ijRr + HVV'llo.Jf), 

k—1 k= 1 

where 0 := max 0^ < 1/2. This finishes the proof. □ 

l^k^d 




max 9 k 
Kl^k^d 
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2.2.2 Property of the nodal interpolation 

Lemma 2.4. For any u £ Pz(K) and v £ Pm{K), it holds that 

A 

hihj 


(V(it - II/ru), Vv)k = ~E E 


d 3 u dv v-^ v-^ 2 3 

- d *+E E 7* h i h * 


3 J K dxidx 2 dx L 

i=l 0 *=1 


I" d 3 u d 3 v 
I K dxidx 2 3 dx 3 


Proof. In order to simplify the calculations, we use auxiliary length = ' , ’" c . 1 ^ i ^ < 

necessary. Since u £ P 3 (K), Taylor expansion yields 

d 

U - " 


tt ft7 <9^1/ / .2* 4 £7 \ 

.-njfU=53 E o + —)• 


2 dxjdxj 

<= 1 * 3 


Thus 


i9(m - IIr-m) _ u u t t- I (>2 4 , r 2\ 


dx; 


J dx 2 dxn 
1 3 


2 dxidx 2 ' 3 3 

l<j<d J 


It follows from the definition of Pm(K ) that 


<9u 


— = Ci,i + C7 2li /ii6 + ^ C' 3 , i (a)( [] /i^r) + 


aeMi 


l^k^d 

k^i 


where Ci,j, C^,*, (^^(a), aridC^,* are all constant coefficients with respect to given u, and Mj is 
multi-indices defined as Mi := {a = (a±, ■ ■ ■ ,ad)\cti = 0 , a*, £ { 0 , 1 }, 1 ^ k ^ i ^ d, |a| > 0 }. 
Elementary calculation yields 

f d^u dv 

J K d^ dx = °> 1 < ^ d ’ 


IK 


hihj d 3 u 
2 dxidx 2 


(£.j -~l+ (c2,ihi£^dx = 0 , 1 < * ^ j ^ d, 


and 


>K 


rfc?-3+ £ ‘)H £ c «<“)( n = 


aeMi 


l^k^d 

k^i 


So, we only need to calculate Ci.,; and C\j, which read 


Ci,i = 


1 


dv 


dx — 


h 2 f d 3 v 


1 d 3 i 


\K\ J K dxi 3|if | J K dx? ^ an< ^ ^ 4,i 2dxP 


Elementary calculation yields 


[ Ci ,i 

JK 


hihj d 3 i 


2 dxidx 2 v 3 3 


fe 2 -o+tf)) dx = - 


C\,ihjhj f d 3 i 


Ik dxidx 2 


dx 


hihj 


d 3 u d 3 v 

i K dx,dx‘j dxC‘ 9 J K dxidx 2 dx 3 


d 3 u dv h 3 hj 
-dx + 13 


dx, 1 ^ i j < d, 


and 


2 t 2 


Ik 


c^ihn 


hihj d 3 i 


2 dxfdxj v 3 3 


(E +c 2 )K = - 


h 3 hi 


d 3 u d 3 v 


15 J K dxfdxj dx 3 
A combination of (12.201) and (12.191) and some elementary calculation yield, 
f d(u — H K u) dv 


dx, 1 ^ i j ^ d. 


/ K dxi dxi 


<* = - J] Mi 


1 <j^d 
j¥* 


IK 


d 3 u dv , v-^ 2 , 3 , 

dx + E 45 ^ 


dxidx 2 dxi 


1 

j¥* 


/iC 


d 3 M d 3 u 
dxidx 2 dx 3 


dx. 


dx. 

i , when 

(2.18) 

(2.19) 

( 2 . 20 ) 

a set of 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 


A summation of (12.261) with respect to i from 1 to d completes the proof. 


□ 
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3 Rectangular Morley element scheme for second order problems 


3.1 Subdivision of the domain and the finite element space 

For simplicity, in this paper, let Q C be a bounded domain which can be subdivided to a rectangular 
triangulation Th- For the triangulation 71, let Mh denote the set of all the vertices, A f h = M' h U M^, with 
Ml and Ml consisting of the interior vertices and the boundary vertices, respectively. Similarly, let 
Th = TI\JT% denote the set of all the (d — l)-faces, with T{ and T\ consisting of the interior faces and 
the boundary faces, respectively. For / 6 T' h , Kj 1 and Kf are the two adjacent elements that share / as 
a common face, and njf and denote the unit outer normal vectors of Kj and Kf, respectively, on /. 
Given a triangulation Th of 12, define the d -rectangular Morley element space 


Vh := {vh £ L 2 (fl) : Vh\x £ Pm(K), Vh(x) is continuous for x £ Mh, / -rr—ds is continuous for / £ Tl}, 

J f dn f 

and associated with 

Vho ■= {vh £ V h : v h (x) = 0, for x £ M^}- 
With respect to the vertices and the faces respectively, we define 

Vh ■■= {Vh £ Vh ■■ f ir^-ds = 0, V/ £ T h }, and Vf := {v h £ V h : v h {x) = 0, Vx £ M h }, 

Jf onf 

and for any f £ Th, 

Vf := {v h £ V h : Vh(x) =0, Vx £ Mh] [ S~^~ds = 0, V/ £ T h other than /}. 

J f' on.fi 


Also, define 


Evidently, 


Vko : = v* n v h0 . 


V h = Vf © V^ = (BferVf © V*, and V h0 = Vf © V* = ® fer vf © V*. 


For each element K £ Th, let hx be the diameter of the smallest ball containing K, and px be the 
diameter of the largest ball contained in K. Let Th belong to a family of triangulations described in 
previous section with h —» 0. We assume that Th satisfied that hx ^ h ^ PPx, VA' £ Th for a positive 
constant r/ independent of h. 

We introduce the following triangulation-dependent norm || ■ \[ m ,h and semi-norm | - \ m ,h'- 


= M 2 m,K) 1/2 A 

KeTh 


i,h = (Y^ m 

K&T h 


,,K 


)l/2 


for such functions v that v\x £ H m (K),VK £ Th- 

3.2 Model problem and its discretization 

We consider the following second order elliptic problem: with / £ L 2 (fl), 

in f 1, 
on dQ, 



(3.1) 


where A is the standard Laplacian operator. 
Define 


du dv 
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The weak form of problem (IT111 is: find u £ Hq (fi) such that 

a(u,v) = (f,v), Vv e 

For v,w £ L 2 (il) that v\k,w\k £ H X {K),MK £ Th, we define 


a h (u,v)= J2 J J2 

KeT h Jk i=1 ° X% ° Xl 


The finite element method for problem (13.11) is: find Uh £ 14o such that 


(3.2) 


ah(uh,v h ) = (f,v h ), Vv h £ V h0 . (3.3) 

Because Vh is continuous on A 4, the weak continuity property ensures the uniqueness of the solution. 

3.3 Some intrinsic properties of 14 

3.3.1 Stable decomposition with respect to vertices and faces 

Firstly, we show that the decomposition of Vh with respect to vertices and faces is stable. 

Lemma 3.1. For any Vh £ 14, there exist uniquely v* £ Vf. vf £ Vf, such that 

v h = v* + vf, 

and moreover, 

\Vh\l,h + \vf\i ,h < \Vh\l,h- 

Proof. Given Vh £ 14, the existence and uniqueness of v* and vf is evident. Now we prove the stability 
of the decomposition. On every cell K , 

(VhVh,VhVh)K = (Y« + Vh),V h (vh +v£))k 

= (y h v% ,V h Vh)K + (V h v £, S7 h vf) K +2(V h v*, W h Vh)K- 

By Lemma T2.31 

l(V/iU*, V h vf) K \ 4 0((V h Vh ,V h Vh)K + {V h vf ,V h vf) K ), 

therefore, 

(VhVh-. V hVh)x ^ (1 — 20)((V/jr>^, VhVh)x + hvf , V/jU^)#). 

Making a summation on every cell K, we obtain that 

(VhV h ,V h v h ) >{l- 2 8)((V h v£,V h v£) + (V h vf,V h v£)). 

Then we have 

l|V h « h ||o,n 0)(\\V h vZ\\o,n + ||V h t^|| 0l n). 

This finishes the proof. □ 


3.3.2 The essential continuity of the finite element functions 

In this section, we estimate the consistency error with respect to different finite element functions. 
Lemma 3.2. The estimate below holds uniformly for any shape-regular family of triangulations. 

1. It holds for v* £ V^ that 

I a h (.v,vf) + (Av,v£)\ < h k ~ 1 \v\k,n\Vh \i,h, Vu £ k = 2,3; 


(3.4) 
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2. It holds for vf G Vf that 


\a h {v,vf) + (Av,vf)\ < h\v\ 2 ,n\vf\ lth , Vv G H 2 {Q). 


(3.5) 


Proof. The assertions fall into the standard finite element analysis. On every element K , let P% : 
L 2 (K) —» P 0 {K), Pf : L 2 (K) —>• Pi{K) be the L 2 orthogonal projection. For the first item, we have 


\ah{v,v%) + (Av,v%)\ 
(by Lemma [2TTTl 


E E [ ^r~ v h n ids\ = \ J2 J2 I - u K v h)mds\ 

jftn t! J sk dxi jr,i Jok dxi 


EE/ 

KeTh i= 1 J 9K dx * 


< 


E E 11^7 - ^ 2 §^\\o,OK\\v^ \\o,OK < h 

KGTh 


K\v\k,K\ V hU,K 


KeTh »=1 
< h k ~ 1 \v\k,n\Vh k = 2,3. 

By the similar technique, for any (d — l)-face /, making use of Lemma [2711 - (2b we obtain (13.51) . The proof 
is completed. □ 

Lemma 3.3. It holds for all such / that ujf is uniform that 

\au f {v,v{) + {Av,v f h )u, f \ < h 2 Uf \v\ 3 , Uf \v f h \i, Uf , Vrfe ^(w/), v G H 3 (uj f ). (3.6) 

Here a Uf (v,v{) = {Vv,S7v{) Kl + (Vv,Vv f h ) KR . 

Proof. By Lemma 12.21 the lemma follows by the same technique as that for Lemma 13.21 □ 

Lemma 3.4. For any d-rectangle K, and / C dK a (d — l)-face, it holds that 

\(Vv,Vv{) K + {Av,v[) k \ < h 2 K \v\ 3 ,k:\vI\i, K , V v f h G Pf(K), v G H 3 {K), v\ f = 0. (3.7) 

Proof. Without loss of generality, assume the normal direction of / to be aq, and / = F\. From the 

f 

expression of the basis function, v J h = 0 on F\ and F<i- Using formula of integration by parts, we have 
| (X7v,Vv f h ) K + (Av,v f h ) K \ = \J2 [ p-v[mds\ 

i=1 JdK ax i 

= 'S ( E ^ Lds ~ L S : vLds)1 = 1 g L U vidA 

Because v\f = 0, = 0(2 < k < d) on /. Then by Poincare inequality, we have ||^|o ,k ^ hi(\v\ 3 ,K- 

Because = 0, 

d 2 i 


IE/ ^A v hdx I = IE/ 

i—2 Jk dX i i =2 Jk dX i 


< 


E Iloilo,- n^|| 0>if < |u| 3 ,if |^|i,J<r- 


This finishes the proof. □ 

Lemma 3.5. The estimate below holds for any shape-regular family of uniform triangulations: 


\a h (v,vf) + (Au,vf)| < h 2 \v\ 3 , Q \v£ Vv G id 3 (U) n ^(H), v h eVf. 


(3.8) 
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Proof. Given vf £ Vf, a decomposition follows that v h 


v ^ h , with v £ Vf. By Lemma 13.31 and 

fer h 


Lemma 13.41 we obtain 

\a h (v,vf) + (Av,vf)\ 



E 1 a h' 

{v,v f h ) + {Av,v[)\ 







< 

<* 2 (E 

\v\3,u f \vl\u f + 

E 

\v\ 3 ,K f \v f h \ Kf ) 


fen 




< 

'> 2 <£ 

\ v \i,uj f + e 

v \i,Kj 

y /2 (E ki?. fc ) 1/a 


ten 





Here, Kf is the element containing (d— l)-face / for / C dfl. Evidently, (X)/ e jr« M§,^ / +S/e. 7 r £ M \,K f ) 1 ^ 2 ~ 
|^| 3 ,n- Besides, by Lemma 1231 


\ v h\i,h = E E ii v ^iio,x< 

E ii v E v h\\o,K = E \\v<\\Ik. 

(3.9) 

f KETh fELTh 

KeTh fGFh K&Th 


Therefore, (13.81) follows. This finishes the proof. 

Lemma 3.6. Given Vh £ 14, there exist uniquely 

v h e V h and v h e V h (/ e Y)> so that 

□ 

V h = v% 

+ E <> 

(3.10) 




and moreover, 



\ v h\i.h + E 

\v f h \l ,h £ \Vh\lh- 

(3.11) 


S^h 

Proof. The proof follows from Lemma 13.11 and (liOl) . □ 

This lemma can hint an optimal fast auxiliary space preconditioner for the discretized system. 


3.4 Convergence Analysis: error estimate in energy norm 

To estimate the convergence rate of the discretization, we begin with the famous Strang lemma below. 

Lemma 3.7. (the second Strang Lemma) Let u and Uh be the solutions of problems (13.11) and (13.21) 
respectively. Then 


\u-u h \i,h= inf \u-v h \i,h+ sup 

Vh&Vh 0 O^WheVho 


\{f, W h ) - a h (u,w h )\ 


\Wh l,h 


Lemma 3.8. Let 14 and 14o be the finite element spaces of the d-rectangular Morley element. Then, 
for k = 2, 3, we have 

k 

inf Y h m \v - v h \ m h < h k \v \ k , Vu £ H k (£l), 


Vh£V h 


m —0 


inf Y h m \v - v h \ m , h < h k \v | fc , Vu £ H k (Sl) n H l Q (Sl). 

vh&Vho ' 
m =0 

Theorem 3.9. Let u and Uh be the solutions of problems (13.11) and (13.21) respectively. Then 

1. if u £ Hq(Q) DiL 2 (H), then 

|u - u h \i,h < h\u\ 2 . 

2. if the triangulation is uniform, and u £ Hq (12) D 7L 3 (H), then 

|u - Uh\i,h < h 2 \u\ 3 . 


(3-12) 


(3.13) 
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Figure 2. Illustration of a general shape regular triangulations in two dimensions. The triangulation in right is a combi¬ 
nation of small patterns as the left one. 


Figure 3 Illustration of uniform triangulations in two dimensions. 


Proof. By Strang lemma and the approximation estimate Lemma 13.81 we only have to study the con¬ 
sistency error, which can be estimated by Lemma |3.21 Lemma 13.51 and the stability of decomposition of 
Lemma HTTl The proof is finished. □ 

Theorem l3.9l reveals that generally on a shape-regular family of triangulations, which can be as “bad” as 
ones shown in Figure[2j the error decays with 0{h) order in energy norm, and on a family of triangulation, 
like one shown in Figure [31 an 0(h 2 ) order can be expected. There is also an intervenient result on 
divisionally uniform triangulation. For a family of divisionally uniform triangulations, we refer to a 
family of conforming triangulations on f l which are uniform triangulations on f tj, j = 1,2,..., J, with 
{£lj}j=i:j a subdivision of SI. (Figured) 


Figure 4. Illustration of divisionally uniform triangulations in two dimension. As shown in the left figure, the domain 
consists of four subdomains, and restricted in each subdomain, the triangulation is uniform. 


Theorem 3.10. Let u and Uh be the solutions of problems (13.11) and (13.21) respectively. If the triangu¬ 
lation is divisionally uniform, and u £ H^fl) D fL 25 (fl), 

\u - Uh\i,h < (3-14) 

For simplicity, we only prove in detail the case that SI is divided to two subdomains. The more 
complicated cases are on the same line. 
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Let 12 = Sli U O 2 , and T = 912 1 0 <9122 be the interface between the two subdomains. A stripe along 
the interface is denoted by 

Qg = {x £ 12|dist(x, T) < 5}. 

The lemma below is important in the technical analysis. 

Lemma 3.11. [2] For v £ where 0 ^ s ^ 0.5, we have: 


IMIoa ;$ <5 S IM| S . 

Lemma 3.12. On divisionally uniform triangulations, it holds for vf £ Vf that 

I ah(v,vf) + (Av,vf) n \ < /i 1,5 ||v|| 2 .s|^ \i,h, v £ H 2 5 {fl) 0 H^(Q). 
Proof. Given vf £ Vf. it can be decomposed to vh = v lr Thus 

a h (v,vf) + (Av,vf) = l + E + E 

\/cnT\r /cfil\r fcr 



It can be proved that 


(3.15) 


E 

/cfii\r 


ah(v,v{) + (Av,v f h ) 


<h k Vlfe.nJI v h\\hh, k = 2,3, 

/ cnT\r 


(3.16) 


E a h( v > v h) + ( Av ’ V h) 

f cfilyr 


<h k 1 |u|fc,n 2 1| J2 v h\\hh, k = 2,3, 
fcTh\r 


and 


E [ a h( V ’ V h) + ( Av , V h) 

fcr 


<h\v\ 2}{ j fcrUlf \\ V h\U,h- 

fcTh\r 


(3.17) 


(3.18) 


Evidently, U/ C rw/ C luh, where H is the size of the biggest cell in and H < h. Therefore, by Lemma 

EU 


l U |2,U 


fc.r 


< h°' 5 \ 


V\\2.5- 


Combining (13.161) . (13.171) and (13.181) leads to (13.151) . This finishes the proof. 


□ 


Proof of Theorem 13.101 The proof follows just the same line as the proof of Theorem 13.91 provided 
Lemma 13.121 □ 

3.5 Convergence analysis: error estimate in L 2 norm 

By standard duality argument, we can prove the upper bound of the L 2 norm of error on convex 
domains. 

Theorem 3.13. Assume 11 is convex, let / £ L 2 (12), and let u and Uh be the solutions of problems (13.11) 
and ed, respectively. Then 

\\u-u h \\ 0 <h 2 \\f\\ 0 . (3.19) 

In general, the estimate above can not be improved. Assuming triangulationes are uniform, we present 
in detail the lower bound estimate of the error in L 2 norm. It follows the idea of m- The main result 
of this section is the theorem below. 

Theorem 3.14. Let u and ug. be solutions of problem ED and ED , respectively. Suppose that 
u £ ifg(12) P|1L S (12), s Sj 3 and s > | + 1. Then, provided ||/||o 7 ^ 0, 

\\u - u h \\ 0 ^/3h 2 , (3.20) 

where /? = 5 is a positive constant, which is independent of the triangulation size h and the 

triangulation size is small enough. 
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Remark 3.15. By the embedding theorem of the Sobolev space, we need higher regularity of the solution 
in higher dimensions in order to guarantee H S (K) C C ll ( K). Furthermore, it ensures the continuity of 
interpolation operators. 

Remark 3.16. For the rectangular domain 12, the condition ||/||o ^ 0 implies that || a ®dx- ||o 7 ^ 0, 

1 ^ * 7 ^ j ^ d. In fact, if \\ q§^ . ||o = 0, 1 < i ^ j < d, then u is of the form u = J2i=i f( x i), for some 
function f(xi ) with respect to Xi. Then the boundary condition indicates u = 0. 

We postpone the proof of Theorem 13.141 after several technical lemmas. 

Firstly, define the global interpolation operator 11/j and P k to V/, by 

^h\K = ^K for any K G Th■ (3-21) 

and 

P k \ K = P k K for any K G T h . (3.22) 

By means of Lemma 12.41 we can obtain the following crucial result. 

Lemma 3.17. For u G Hq (12) f) R s (12), s ^ 3 and s > | + 1, it holds that, 

(V h {u - U h u), Vh^hu) ^ ah 2 , (3.23) 


for some positive constant a, which is independent of the triangulation size h provided that ||/||o 7 ^ 0 
and that the triangulation size is small enough. 

Proof. Given any element K, we follow the idea of m to define Pkv G Ps{K) by 


/ V l Pxvdx= / X7 l vdx, l = 0,1,2,3, 
Jk Jk 


(3.24) 


for any v G H s ( 12), (s ^ 3 and s > | + 1). Note that the operator Pk is well-defined. The interpolation 
operator Pk has the following error estimates: 


|u — Pkv\j,k ^ h 3 i\v\ 3 t K, j = 0,1, 2, 3, 
|n - Pkv\j,k < h\v\j + i,K, j = 0,1,2, 


(3.25) 


provided that v G LP(12), (s > 3 and s > | + 1). It follows from the definition of Pk in (13.241) that 


\7 3 P k v = P&VV 

By the aid of Pk, we have the following decomposition 

(V/ t (zi — Il^n), V/jIIfcu) = ^ (V/i(Pjcm - HkPku), V/iII^m)#: 


(3.26) 


KGT h 


+ ^2 (Vh(Id — H K )(Id — Pk)u, \/hJd K u)K 
KeT h 

= J\ + J2- 

By means of Lemma [^H the first term Ji on the right-hand side of (13.271) can be rewritten as 


(3.27) 


* = -£££ 


KeTh i= 1 1 


hihj f 8 3 Pku OHku 
dxidx] dxi 


IK 


di +EE E 45^. 


KeTh i= 1 

3+i 


IK 


d 3 P K u d 3 U K u 
dxidx 2 dx 3 


da: 


= -£££ 


hihn 


d 3 u du 


d 


dx + E E E 


KsThi= ll^d 3 JKdx^dXi - ■ 3 


hihj j' d 3 {Id — Pk)u OTIku 
dxidx 2 dxi 


da; 


ik 


EE E 

KGTh i=l 1 


hihj f d 3 u d(Id — Hk ) 


d 


l K dxidx j dxi 


di + E E E 45^^ 


KeTh i=l 


r d 3 P K ud 3 n K u 

Ik dxidx 2 j dx 3 


dx. 
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Since the triangulation is uniform and the boundary condition u = 0, on 9f2, thus -§jr\r x = ( 
r Xj . is the face of dQ perpendicular to Xj axe and j ^ i, integrating by parts yields 


EEEt 

KeTh *=1 1 


f d 3 u du 
I K dxidx'j dxi 


da; 


_ Sy^ hihj 


KeTh i=l l^j^d 


f / d 2 u >.2 . h^hj f d 2 u du 

h dxidxj X 2—/ ^ 3 J T dxidxjdxi 

K 2 ffeT/, i=i lx i 2 


ds 


= - EE E 

KEl~h i=l l^j^d 


hihj 11 9 2 u 11 2 

3 ^ dx^Xj ^ 0 ’- K ” 


By the commuting property of (13.261) . 

d 3 {Id — P K )u 


dxidxj 


= (Id-P° K ) 


d 3 i 


dxidx 2 


, 1 < * 7 ^ j < d- 


Note that 


E 


I C®if U I I 

I ar, llo.-K' 


< M 3 ,if- 


This and the error estimate of (13.251) yield 


*=EEE i j i lla^-lll K + o(* 2 )ll(M-n)vM| 0 , K | u |3, K . 

KeTh i=i i<i<d 1 2 

We turn to the second term J 2 on the right-hand side of (13.271) . By the Poincare inequality, 
commuting property of (13.261) . 


\j* I < h 2 £ ||V^(/d-P Jf ) U || 0 , Jf H3</l 2 ||(/d-P 0 MT i || 0 H3. 
KdTh 

Since the piecewise constant functions are dense in the space L 2 (fl), 

\\(Id — P°)V 3 u||o — > 0, when h —> 0. 

Summation of (13.281) . (13.291) and (13.301) completes the proof. 

Again, the lemma below can be found in m- 
Lemma 3.18. Let u and Uh be solutions of problem (13.11) and (13.21) . respectively. Then, 


(-/, u-u h ) = a h (u, U h u - u h ) - (/, - u h ) 

+ a h (u- II h u, u - n h u) + a h {u- II^u, u h - rh.it) 
+ 2(/, II h u -u) + 2 a h {u - U h u, rtft.it). 


Proof of Theorem 13.141 It follows from Lemma l3T2l that 

a h (u,U h u - u h ) - (/, Ib.it - u h ) < h 2 \u\ 3 \U h u - u h \ h < /i 4 |u|g. 
By the Cauchy-Schwarz inequality and the error estimate Theorem 13.81 it yields 
a h {u - U h u,u-Il h u) + 2(f,U h u - u) < h 4 (\u\ 3 + ||/|| 0 )|u| 3 , 
a h (u - II h u, u h - II h u) < h A \u\l- 


, where 


(3.28) 

and the 

(3.29) 

(3.30) 

□ 

(3.31) 

(3.32) 

(3.33) 

(3.34) 































16 


Meng, Yang and Zhang Sci China Math for Review 


The error estimate of the last term of (13.3111 by Lemma 13.171 gives 

ah 2 ^ ah(u — II h u, II hu). (3.35) 

Hence, a combination of (13.3111 - (13.351) leads to 

(-f,u-u h )^Sh 2 . (3.36) 


for some positive constant <5, which is independent of the triangulation size h and the triangulation size 
is small enough. 

Therefore, 


\u-Uh\\o= sup 

o^weL 2 (n) 


(w, u-u h ) > (-/,u- u h ) > 


Iloilo 


-/ Ik 


-h 2 . 


(3.37) 


This finishes the proof. 


□ 


4 Numerical example 

In this section, we present some numerical results of the two-dimensional and three-dimensional RM 
element by uniform triangulation, divisionally uniform triangulation, and general shape regular triangu¬ 
lation, respectively, of domain f l to demonstrate our theoretical results. We follow the approaches shown 
in Figures O [3] and [2] to generate divisionally uniform triangulations, uniform triangulations, and general 
shape regular triangulations, respectively. 

4.1 Two-dimensional examples 

For two-dimensional experiments, we choose the computation domain to be H = [0, l] 2 . We choose / 
such that the exact solution is u\(x,y ) = x(l — x)y( 1 — y) and U 2 {x,y) = sin(7ra;)sin(7r?/), respectively. 
We run the numerical experiments with respect to different kinds of triangulations, and record the 
convergence rate in Figures [5] (for uniform triangulations), [6] (for divisionally uniform triangulations), 
and [7] (for general shape regular triangulations), respectively. 


The errors of u^x.y) The errors of u 2 (x,y) 




Figure 5 The errors in L 2 and H 1 norms with respect to ui(x,y) and U 2 (x,y) on uniform triangulations. 
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The errors of u 1 (x,y) The errors of u 2 (x,y) 



Figure 6 The errors in L 2 and H 1 norms with respect to 

The errors of u 1 (x,y) 



□ HI err 
O L2err 




-8 - 1 - 1 - 

12 3 4 

log(DOFs) 


u±(x,y) and U 2 (x,y) on divisionally uniform triangulations. 

The errors of u 2 (x,y) 



Figure 7 The errors in L 2 and H 1 norms with respect to ui(x,y) and U 2 (x,y ) on general shape regular triangulations. 


4.2 Three-dimensional examples 

For three-dimensional experiments, we choose the computation domain to be In Q = [0, l]' 5 . We 
choose / such that the exact solution is u\(x,y,z) = x{l — x)y{\ — y)z{ 1 — z) and U 2 (x,y,z) = 
sin(7ra;)sin(7r?/)sin(7rz), respectively. We run the numerical experiments with respect to different kinds of 
triangulations, and record the convergence rate in Figures [8] (for uniform triangulations), |9] (for division¬ 
ally uniform triangulations), and [10] (for general shape regular triangulations), respectively. 

5 Concluding remarks 

In this paper, we present a rigorous analysis of the RM element applied for second order problem in 
arbitrary dimensions. To be combined with the standard framework, some special properties of the RM 




































18 


Meng, Yang and Zhang Sci China Math for Review 


The errors of u 1 (x,y,z) The errors of u 2 (x,y,z) 



Figure 8 The errors in L 2 and H 1 norms with respect to m(x,y,z) and U 2 (x, y, z) on uniform triangulations. 


The errors of u 1 (x,y,z) The errors of u 2 (x,y,z) 



Figure 9. The errors in L 2 and H 1 norms with respect to u\{x , y , z) and U 2 (x , y , z) on divisionally uniform triangulations. 


element functions are revealed and used. Both the energy norm and the L 2 norm of the error are studied, 
and the upper bound and the lower bound obtained illustrate that the analysis presented here is optimal. 
The RM element pretends to be one fit for fourth order elliptic perturbation problems, and can also 
expect application for contact/obstacle problems (c.f., e.g., [TIEISI ) in the future. 

The fundamental role of stable decomposition for implementing the “divide and conquer” strategy is 
corroborated again in the analysis of the error estimate in energy norm. Also, the stable decomposition 
as revealed by Lemma 13.61 can be used to design an optimal preconditioner under the framework of fast 
auxiliary space preconditioning |13II46U471I54| . Besides, it is quite interesting to note that the RM element 
space does not contain a nontrivial conforming subspace. This unusual fact makes the a posteriori error 
analysis of the scheme a problem which absorbs theoretical interests. These will be studied in future 
works. 
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Figure 10. The errors in L 2 and H 1 norms with respect to u±(x, y , z) and U 2 (x, y , z) on general shape regular triangulations. 
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